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We investigate the properties of a homogeneous dipolar Bose gas in a weak three-dimensional 
(3D) isotropic speckle disorder at finite temperature. By using the Bogoliubov theory (beyond 
the mean field), we calculate the condensate and the superfluid fractions as a function of density 
and strengths of disorder and interaction. The influence of disorder on the anomalous density, the 
chemical potential and the ground state energy is also analyzed. We show that the peculiar interplay 
of the DDI and weak disorder makes the superfluid fraction and sound velocity anisotropic. 
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I. INTRODUCTION 


Disordered Bose gas in a weak random external po¬ 
tential (dirty boson) represents an interesting model for 
studying the relation between Bose-Einstein condensa¬ 
tion (BEC) and superfluidity and has been the subject 
of many theoretical investigations in the last two decades 
M- ' Experimentally, the dirty boson problem, was 
first studied with superfluid helium in aerosol glasses 
(Vvcor)fTol4T2j. Recently, several groups have loaded ul¬ 
tracold atoms into optical potentials and studied BECs 
in the presence of disorder [14l - [l9l |. 


What happens to a homogeneous BEC if a weak ran¬ 
dom external potential is switched on? Indeed, the pres¬ 
ence of disordered potential may lead to decrease both 
BEC and superfluidity. Furthermore, one of the intrigu¬ 
ing feature of disordered Bose gas is the appearance 
of the so-called Anderson localization [2l| in the non¬ 
interacting case. This phenomenon which can be un¬ 
derstood as the effect of multiple reflections of a plane 
wave by random scatterers or random potential barriers, 
has recently attracted a great deal of interest 
Experimentally, the random potential can be created us¬ 
ing different techniques, one of which is the static laser 
speckle, whereas potential felt by atoms is proportional 
to the speckle intensity with the sign of the detuning 
from the atomic transition 22], Laser speckle, produced 
by passing expanded laser beam through diffusive plates, 
are special in that they have (i) exponential, i.e. strongly 
non-Gaussian, intensity distribution and (ii) finite sup¬ 
port of their power spectrum [23|. Recent progress in 
different experimental realizations of laser speckle disor¬ 
der is reported in references [22|, [24[ . 

In their recent work Abdulaev etal [25}, have shown 
that a Gaussian approximation of the autocorrelation 
function of laser speckles, used in some recent papers, is 
inconsistent with the general background of laser speckle 
theory. They also pointed out that the concept of a quasi- 
3D speckle, which appears due to an extension of the au¬ 
tocorrelation function in the longitudinal direction of a 
transverse 2D speckle, is not applicable for the true 3D 
speckle, since it requires an additional space dimension. 
In this context, they derived an appropriate autocorrela¬ 


tion function for an isotropic 3D laser speckle potential 
which has the Fourier transform given in Ea. dl^l) (see 
below). 

Recent progress in the physics of ultra-cold gases have 
led to the creation of BECs with dipole-dipole interaction 
(DDI ) and stimulated a tremendous boost in theoreti¬ 
cal and experimental studies of weakly interacting Bose 
gases j26l-l28l|. What is important in such systems is that 
the atoms interact via a DDI that is both long ranged 
and anisotropic. By virtue of this interaction, these sys¬ 
tems are expected to open fascinating prospects for the 
observation of novel quantum phases in ultracold atomic 
gases. On the other hand, dipolar BECs confined in ran¬ 
dom media remain largely unexplored. One can quote 
for example, uniform dipolar Bose gas with a Gaussian 
disorder correlation function, a Lorentzian, and a delta- 
correlated disorder have been explored recently by Pel- 
ster etal [29l - l3l| . 

In the present paper, we study the impact of a weak 
disorder potential with 3D isotropic laser speckle auto¬ 
correlation function of Ref. [25] on the properties of a ho¬ 
mogeneous dipolar Bose gas at finite temperature. To 
this end, we use the Bogoliubov theory (beyond the mean 
field) and we calculate in particular the condensed deple¬ 
tion and the anomalous fraction. This latter quantity 
which grows with increasing interactions and vanishes in 
noninteracting systems [32H35I ]. is important to fully un¬ 
derstand the interplay of disorder and interactions. We 
show, in addition, how the anisotropy of the DDI enhance 
quantum, thermal and disorder fluctuations as well as the 
superfluid fraction. 

The rest of the paper is organized as follows. In SecITTl 
we describe our model of the dipolar dilute Bose gas in 
a general random potential. In Sec lIIIl we derive ana¬ 
lytical expressions for the condensate fluctuations and 
some thermodynamic quantities for 3D isotropic laser 
speckle disorder potential at finite temperature. We show 
that the competition between both contact interaction- 
disorder and DDI-disorder leads to enhance the conden¬ 
sate depletion, the anomalous density, disorder fluctua¬ 
tion, ground state energy, equation of state and the sound 
velocity. In Sec lIVl the superfluid fraction is obtained 
and its characteristics are discussed. Finally, our conclu- 
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sions and outlook remain in SecfVl 


II. THE MODEL 


where b k and 6k are operators of elementary excita¬ 
tions. The Bogoliubov functions u k ,v k are expressed in 
a standard way: u k ,v k = (y/ekjEk ± \JE k /e k )/2 with 
E k = h 2 k 2 /2m is the energy of a free particle, and 


We consider the effects of an external random field 
on a dilute 3D dipolar Bose gas with dipoles oriented 
perpendicularly to the plane. The Hamiltonian of the 
system is written as: 


^ = M^rU k . (6) 

V v 

The Bogoliubov excitations energy is given by 


H = y d 3 r ^-^—A + C/(r)^ ip(r) 

+ iy d 3 r J d 3 r'ft (r)ft (r')V(r — r')ip(r')ip(r), (1) 


where ij}' and if) denote, respectevly the usual creation 
and annihilation field operators, the interaction potential 
V(r — r') = g8(r — r') + Vd d ( r — r'), g = 4nfi 2 a/m corre¬ 
sponds to the short-range part of the interaction and is 
parametrized by the scattering length as a. On the other 
hand, the dipole-dipole component reads 


V d {r) 


C d d 1-3 cos 2 9 
47t r 3 


( 2 ) 


where the coupling constant Cdd is MqM 2 for particles 
having a permanent magnetic dipole moment M (Mo is 
the magnetic permeability in vacuum) and d 2 /eo for par¬ 
ticles having a permanent electric dipole d (eo is the per¬ 
mittivity of vacuum), m is the particle mass, and 9 is 
the angle between the relative position of the particles r 
and the direction of the dipole. The characteristic dipole- 
dipole distance can be defined as r* = mCdd/^nh 2 . For 
most polar molecules r* ranges from 10 to 10 4 A. The 
disorder potential is described by vanishing ensemble av¬ 
erages (U( r)} = 0 and a finite correlation of the form 
(U (r)U(r')) = R{ r,r'). 

Passing to the Fourier transform and working in the 
momentum space, the Hamiltonian (pQ) takes the form: 


a r—> Ti k^ + ^ 1 ^ A , 

H = / v 2 ^ a k Q k + ~p/ Uk-p a k a p (3) 

k k,p 

+ 2Y ^ (P)“k+q“k-q“ k -T“ k -P’ 
k.q.p 


where V is a quantization volume, and the interaction 
potential in momentum space is given by [35| 

/( k ) =ff[ 1 + e dd (3cos 2 6»-l)], (4) 


Assuming the weakly interacting regime where r* -C £ 
with £ = h/y/mgn being the healing length and n is the 
total density, we may use the Bogoliubov approach. Ap¬ 
plying the inhomogeneous Bogoliubov transformations 

jf: 

dk = Mfc6k-Dc6l k -/3k, a k = u fc 6j i -u fc 6_k-/3k, (5) 


e k = y/E 2 + 2^ od {0)E k , (7) 

where fiQd = n lim /(k) is the zeroth order chemical po- 
fc—»o 

tential. 

Importantaly, the spectrum 0 is independent of the ran¬ 
dom potential. This independence holds in fact only in 
zeroth order in perturbation theory; conversely, higher 
order calculations render the spectrum dependent on 
the random potential due to the contribution of the 
anomalous terms (see below). For k —> 0, the excita¬ 
tions are sound waves e k = hc s d(9)k , where c s d{9 ) = 
c s y/l + edd (3 cos 2 6—1) with c s = \Jgn/m is the sound 
velocity without DDL Due to the anisotropy of the dipo¬ 
lar interaction, the sound velocity acquires a dependence 
on the propagation direction, which is fixed by the angle 
9 between the propagation direction and the dipolar ori¬ 
entation. This angular dependence of the sound velocity 
has been confirmed experimentally [36j . 

Therefore, the diagonal form of the Hamiltonian of the 
dirty dipolar Bose gas m can be written as 

H = E+Y / £kblb^ (8) 

k 


where E = Eod + dE + Er : 

Eod(9) = g,od(9)N/2 with N being the total number of 
particles. 

«?=i^[e fe -£ fe -n/(k)], (9) 

k 

is the ground-state energy correction due to qunatum 
fluctuations. 


E R = -Y j n{\U k \ 2 )^ = ^Y.nRk^t, (10) 

gives the correction to the ground-state energy due to the 
external random potential. 

The noncondensed and the anomalous densities are de¬ 
fined as h = J]k(®k®k) and m = )T) k (dkd-k)j respec¬ 
tively. Then invoking for the operators a k the trans¬ 
formation ©, setting (6 k 6 k ) = d^^Nk and putting 
the rest of the expectation values equal to zero, where 
N k = [exp (e k /T) — l]^ 1 are occupation numbers for 
the excitations. As we work in the thermodynamic 
limit, the sum over k can be replaced by the integral 
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= V f d 3 k/{2n) 3 and using the fact that 2N(x) + l = 
coth(a:/2), we obtain: 


n = — 
2 


r d 3 k 

E k + /(k )n 

coth ( £k ^ 1 

j (2tt) 3 

£k 

1. 1 V2T/ 1 \ 


+n R , (11) 


and 


in = — 


1 

2 


d 3 k /(k)n 
(2tt)3 


coth 



+ n R . 


( 12 ) 


The contribution of the random potential comes through 
the last terms in Eqs m and m- These terms are 
dehned as 


n R = = n 

k 


d 3 k E\ 

( 2^) 3 4 


Rk- 


(13) 


Expressions CD and CD must satisfy the equality 
** (Rl + 11 “ '“‘ |2 = J s inh i ( £t /2T) 

(14) 

Equation (fT4l) clearly shows that fh is larger than n at 
low temperature irrespective to the presence of an ex¬ 
ternal random potential or not. So the omission of the 
anomalous density in this situation is principally unjus¬ 
tified approximation and wrong from the mathematical 
point of view }33l435j . 


III. BEC FLUCTUATIONS AND 
THERMODYNAMICS QUANTITIES IN WEAK 
ISOTROPIC LASER SPECKLE 

To proceed further in practical calculations, we must 
define the laser speckle potential: U(r) = Uo + AU(r), 
where Uo is defined by the light far-held intensity as Uo = 
(I) and (AU(r)) = 0. At the derivation of U{ r), it was 
assumed that the incident laser wave does not induce an 
atomic electron interlevel transition, but merely deforms 
the atomic ground state. It is useful now to specify the re¬ 
lationship between the far-held intensity autocorrelation 
function |C/(r)| 2 , the laser speckle autocorrelation func¬ 
tion |CU(r)| 2 and the disorder potential correlation func¬ 
tion. One can write then: |C/(r)| 2 = (U(r')U(i J + r)} 
and \Ca{4\ 2 = {AU(r')AU(r' + t))/Uq. Therefore, us¬ 
ing the Fourier transform, we get 

|CV(k)| 2 = [7g [<5(k) + |CU(k)| 2 ], (15) 

where the autocorrelation function of the laser speckle is 
given by [25| 

|CU(k)| 2 = -^-(2tr) 3 [(2crfc) 3 - 12(2ak) + 16], (16) 


where cr characterize the correlation length of the dis¬ 
order (for further computational details, see Ref |25j). 
Interestingly, we see from the formula of |C( 4 (k)| 2 that 
its value becomes zero for k = 1/a. Hence, the momen¬ 
tum in (1161) only varies in a hnite interval from zero, in 
contrast to the case for a Gaussian function^. 

Putting l?(k) = l?|CU(k)| 2 [25}], where R=Uq stands 
for the disorder strength. Substituting now the function 
(1161) in equation (1131) and performing the integration over 
the momentum form 0 to 1 /cr, we get the expression for 
the condensate fluctuation due to the external random 
potential 


n R 


m 2 R 

87T 3 / 2 ?! 4 


h(edd, g:), 


where 

h(e dd , a) = [ dO smOS^al 

Jo y/ 1 + £ dd( 3cos 2 6-1) 

is depicted in Fig ID and the function 


(17) 


(18) 


rt ( \ 1 /D" 

S(a) = 2iV2 


4 — (8a + 6) In 1 H- 

2a 


arctan 


(,V2a) 


with a = cr 2 [l + £dd (3 cos 2 9 — l)]/£ 2 . 



FIG. 1. (Color online) Behavior of the disorder function 
h(tdd, c/O from equation (1181) . as a function of <r/£. Red 
line: Er atoms (tdd, = 0.38). Blue line: pure contact interac¬ 
tion (tdd = 0). 

In the absence of the DDI (edd = 0), we recover the 
result for the 3D BEC with short-range interparticle in¬ 
teraction of Ref [25|. For cr/£ —> 0 and e dd = 0, we 
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read off from Eg. flTSl) that one obtains h(edd,&) —> 1 
(see also FigHJ). Therefore, we should reproduce the 
Huang and Meng result [If for the condensate depletion 
in this limit. For cr/£ —> 0, we get from Eg. flTSl) that 
h{edd, 0) = Q-i(edd)- Thus, the disorder fluctuation flTTl) 
becomes identical to that obtained in 3D dipolar BEC 
with delta correlated disorder [31 ] 


n _h 


m 2 R 

8ir 3 / 2 h 4 


Q~l(^dd), 


(19) 


where the contribution of the DDI is expressed by the 
functions Qj(e d d) = (1 - £dd) 3/2 ?Fi (- 5 , 5 ; , 

where 2 F 1 is the hypergeometric function. Note that 
functions Qj(e d d) attain their maximal values for edd ~ 1 
and become imaginary for edd > 1 [HSzi. 

On the other hand, the function m decreases with in¬ 
creasing disorder correlation length while it rises for in¬ 
creasing edd and diverges in the limit edd > 1- Another 
important consequence is that when a vanishes, hr be¬ 
comes infinite. This means that the system would col¬ 
lapse if there were no repulsive interactions between par¬ 
ticles. 

Upon calculating integral in Eq. m, we get for the 
condensate depletion 


n 

n 


8 

3 


na 3 / x 2 

- QoK^dd) + o 

7T O 



2 

Q-l(tdd) 


+ 2nR' 


~h(e d d, ^), 


( 20 ) 


where R' = R/g 2 n is a dimensionless disorder strength. 

The integral in Eo. flTSl) is ultraviolet divergent. This 
divergence is well-known to be unphysical, since it is 
caused by the usage of the contact interaction potential. 
A general way of treating such integrals is as follows. 
First, one restricts to asymptotically weak coupling and 
introduces the Bcliaev-type second order coupling con¬ 
stant [35:] 


at higher temperatures i.e. T gn , where the main con¬ 
tribution to integrals m and m comes from the single 
particle excitations. Hence, the thermal contribution of h 
becomes identical to the density of noncondensed atoms 
in an ideal Bose gas [35| . while the thermal contribution 
of rh tends to zero since the gas is completely thermalized 
in this range of temperature [12, [12, HH ■ The last term 
in (12U1) and (1221) accounts for the effect of disorder on the 
noncondensed and the anomalous densities. 

Equation (1221) clearly shows that at zero temperature, 
the anomalous density is larger than the noncondensed 
density for any range of the dipolar interaction as well as 
for any value of the strength and the correlation length 
of the disorder as it has been anticipated above. More¬ 
over, rh changes its sign with increasing temperature in 
agreement with uniform Bose gas with a pure contact 
interaction [35[ . Likewise, the anomalous density ob¬ 
tained in (1221) permits us to elaborate in a straightforward 
manner the equation of state and thus, leads to a finite 
compressibility (see below). Remarkably, Eos. (12Ul) and 
(1221) reproduce the short range interaction results since 
Qj(edd = 0) = 1. Furthermore, the DDI enhances quan¬ 
tum, thermal and disorder fluctuations of the condensate 
for increasing edd as is shown in FigQ] 

The Bogoliubov approach assumes that fluctuations 
should be small. We thus conclude from Eqs. 
and (1221) that at T = 0, the validity of the Bogoli- 
ubov theory requires inequalities Vna 3 Qo(edd) ■C 1 and 
RWna 3 h(e dd , a) <C 1. For R' = 0, this parameter differs 
only by the factor Qo(edd) from the universal small pa¬ 
rameter of the theory, V na 3 <C 1, in the absence of DDI. 
At T <C gn , t he Bogoliubov theory requires the condition 
{T/ gn)Vncfi Q-\{edd) 1. The appearance of the extra 
factor ( T / gn) originates from the thermal fluctuations 
corrections. 

The presence of quantum and disorder fluctuations 
leads also to corrections of the chemical potential which 
are given by Sg = X) k /( k )[ w fc( v fc - u k )] = D k /( k )(« + 
rh) [ill, uMI. Inserting the definitions (Jill) and flT2l) 

into the expression of Sg, we find after integration: 


/*(k) = /( k ) - ^ 


d 3 q /(-q)/(q) 
(2 tt)3 2 E q 


( 21 ) 


After the subtraction of the ultraviolet divergent part, 
the anomalous fraction turns out to be given 


m 

n 


—a(^) -1 

7 r 3 


na° 
7 r 


ttT 

gn 


Q-li^dd) 


+ 2t tR’ 


-h(e d d,u). 


( 22 ) 


The leading term in Eos. (EUl) and (E^l) represents the 
qunatum fluctuation[35|. The subleading term which 
represents the thermal fluctuation[35}, is calculated at 
temperatures T <C gn, where the main contribution to 
integrals m and m comes from the region of small 
momenta (e*, = hc s dk). The situation is quite different 


Sg 

go 


32 

y 


na° 


na° 


Qb(,f-dd) + 4:TvR'\ - hi (edd, Oi). (23) 


where h 1 (e d d,ot) = f* dO sin 0-^/1 + edd(3cos 2 0 - l)S(a) 
and go = gn. 

In the absence of disordered potential (/?' = 0), 

Eq. (1251) coincides with that derived recently in [22, |37( . 
For a condensate with a pure contact interaction 
(Q 5 (edd = 0) = 1) and for R' = 0, the obtained cor¬ 
rection to the chemical potential (1231) excellently agrees 
with the seminal Lee-Huang-Yang quantum corrected 
equation of state @. 

The energy shift due to the interaction and the quan¬ 
tum fluctuations J9| is ultraviolet divergent. The diffi¬ 
culty is overcome if one takes into account the second- 
order correction to the coupling constant ©. a simple 
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calculation yields [35, 37] 

6E=^Vgn 2 J^Q 5 (e dd ). (24) 

However, the energy shift due to the external random 
potential m is evaluated as 

^ = IQ-kR' \E^—hi(e dd ,a), (25) 

-CjQ V 7T 




where Eq = Ngn/ 2. 

When a <C £, the energy shift due to the external ran¬ 
dom potential m is ultraviolet divergent. Again, by 
introducing the renormalized coupling constant m one 
gets: E r /E 0 = IQttR' x/na 3 /nQi(edd) which well con- 
cides with the result obtained with delta correlated dis¬ 
order of Ref [lU . 


IV. SUPERFLUID FRACTION 


The superfluid fraction n s /n can be found from the 
normal fraction n n /n which is determined by the trans¬ 
verse current-current correlator n s /n = 1 — n n /n. We 
apply a Galilean boost with the total momentum of the 
moving system P = mti(nv s + n n v n ), where v s denotes 
the superfluid velocity and and v n = u — v s is normal 
fluid velocity with u being a boost velocity j3Cj. The 
superfluid fraction is then written 


— — S 

- '■'7.' 


-4 


2 

Tn 


d 3 k Ti 2 nRkkikj 
(27t) 3 2m E k [E k - 2?r/(k )] 2 
d 3 k h 2 kikj 
(27t) 3 _2m4sinh 2 (e fc /2T). 


(26) 


It is worth remarking that if in expression (12611 rh were 
omitted, then the related integral would be divergent 
leading to the meaningless value n s —> 00 . This indicates 
that the presence of the anomalous density is crucial for 
the occurrence of the superfluidity in Bose gases @, [3j| 
which is in fact natural since both quantities are caused 
by atomic correlations. 

Equation. (l26l) yields a superfluid density that depends 
on the direction of the superfluid motion with respect to 
the orientation of the dipoles. In the parallel direction, 
the superfluid fraction reads 


ns . . ina 3 u r 2 h 

— = I-471T? \ - h n (e dd ,a)- 

n ~ 


where the function 


hH e dd,a) = / 
Jo 


d9 


45miic. \ he. 


sin 8 cos 2 OS (a) 


) Q-5( e dd), 


0 \/l + edd {3 cos 2 0 - 1 ) ’ 


(27) 

(28) 


is decreasing with increasing edd for fixed cr/f; as is de¬ 
picted in Figj2]a. And the functions Q^{e dd ) = |(1 — 


FIG. 2. (Color online) Behavior of the disorder functions 
h n (e dd ,a/^) (a) and h ± (e dd ,a/^) (b), as a function of er/£. 
Red line: Er atoms (edd = 0.38). Blue line: pure contact 
interaction ( t dd = 0 ). 


edd) j/2 2 F 1 (-§, |; §; have the following proper¬ 

ties: Q^(e dd = 0) = 1/3 and imaginary for e d d > 1 HH- 
In the perpendicular direction, the superfluid fraction 
(l26l) takes the form 


= 1-2tt R‘ 


na 3 1 . . n 2 h 

- h {e dd ,a)-— - 

n AbmnCs 


T 

hc s 


where the function 


Q^sitdd), 

(29) 


h ± (e dd ,a)=[ d0^ 
Jo V 1 


sin 0(1 — cos 2 6)S(a ) 


e dd {3 cos 2 0 - 1 ) 

= h{e dd ,a) - h^(e d d,a), (30) 


is increasing with e dd for fixed cr/£ as is displayed in 

FiglUb. And Qf{e dd ) = Qj(tdd) - Q^(e dd )- 

Third terms in (l27l) and (l29l) which represent the ther¬ 
mal contribution of njr and nl , are calculated at low tem¬ 
peratures T -C ng. Whereas, at T> ng, there is copious 
evidence that both thermal terms of n s concide with the 
noncondensed density of an ideal Bose gas. Furthermore, 
we read off from Eos. (l27l) and (l29l) that for e dd < 0.5, the 
thermal contribution of njr is smaller than that of nl, 
while the situations is inverted for e dd >0.5. 

For <t/£ —>■ 0 and e dd = 0, both components of the 
superfluid fraction reduce to n s /n = 1 — 4n/{/3n, which 
well recoves earlier results of Refs @,il for isotropic 
contact interaction. For a/l; —> 0, we have h^{e dd , 0) = 

Qli(edd) and h ± (e dd , 0) = Q/q^d). As a result, the 
disorder correction to superfluid fraction o becomes 
identical to that obtained in 3D dipolar BEC with delta 
correlated disorder [3l| . 

We should stress also that for increasing e d d , h)\ (e dd , a) 
decreases, whereas h(e dd ,ct) increases for fixed a /£. 
Therefore, this reveals that there exists a critical value 
of interaction e dd beyond which the system has the sur¬ 
prising property that the disorder-induced depletion of 
the parallel superfluid density is smaller than the con¬ 
densate depletion even at T = 0. This can be attributed 
to the fact that the localized particles can not contribute 
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to superfluidity and, hence, form obstacles for the su¬ 
perfluid flow. For large disorder correlation length i.e. 
a £, e dd reduces indicating that the localized particles 
are localized in the respective minima for the disorder 
potential only for a finite localization time [40j. This lo¬ 
calization time remains to be analyzed in more detail in 
a future work. In addition, the superfluid fraction can 
be either larger or smaller than the condensate fraction 
n c /n = 1 — n/n, depending on temperature, on interac¬ 
tion and on the strength of disorder. Increasing R' leads 
to the simultaneous disappearance of the superfluid and 
condensate fractions. 


Note that the sound velocity of a dipolar BEC in a 
weak external disorder potential can be calculated within 
the hydrodynamic approach as c^(q) = ( d/i/mdn)ci T h s q 
[29l [ 30 ]. where the tensorial property of the superfluid 
density has been taken into account. From Eqs. m and 
it follows that the sound velocity can also be sep¬ 
arated into a parallel and a perpendicular components. 
Both components change via effects of the interaction 
strength e dd , disorder strength R' and the ratio cr/£. One 
can easily show also that the sound veclocity is consi- 
tent with the inverse compressibility kV 1 11 = n 2 dn/dn\ 4l|, 
where the increase in k _ 1 tends to increase the sound ve¬ 


locity and vice versa. 

V. CONCLUSION 

In this paper, we have studied the properties of a ho¬ 
mogeneous dipolar Bose gas in the presence of a weak 
disorder with autocorrelation function for an isotropic 3D 
laser speckle potential at finite temperature. Using the 
Bogoliubov approach, we have calculated the condensate 
fluctuation due to disorder, as well as the correspond¬ 
ing corrections the condensed depletion, the anomalous 
fraction, the chemical potential, the ground state energy 
and the sound velocity. We have pointed out that the 
interplay between the anisotropy of DDI and the exter¬ 
nal random potential leads to modify both the BEC and 
the superfluidity characteristics. Furthermore, we have 
reproduced the expression of the condensate fluctuations 
and thermodymanics quantities obtained in the literature 
in the absence of both the DDI and the disordered poten¬ 
tial. We discuss the validity criterion of the Bogoliubov 
approach in a dirty dipolar BEC. 

Finally, an interesting question that begs to be asked 
is how the interplay of disorder and DDI can affect An¬ 
derson localization, or the quantum phases that arise due 
to disorder in the regime of strong correlations. 
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